Relaxation of Hot Quasiparticles in a d-Wave Superconductor 
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Motivated by recent pump-probe experiments we consider the processes by which "hot" quasi- 
particles produced near the antinodes of a d-wave superconductor can relax. We show that in a 
large region of momentum space processes which break Cooper pairs are forbidden by energy and 
momentum conservation. Equilibration then occurs by scattering with thermal quasiparticles: Um- 
klapp scattering is exponentially suppressed at low temperatures, but small-angle scattering leads to 
power-law behavior. By solving the Boltzmann equation analytically we make detailed predictions 
for the temperature and intensity dependence of these processes, which we compare with experiment. 
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Introduction. — Despite intensive effort for more than 
a decade, no consensus exists on the form and origin 
of the effective electron-electron interaction leading to 
pairing in the cuprate superconductors. Furthermore, 
although it is generally agreed that the superconducting 
state has d-wave symmetry, many properties of the Bo- 
goliubov quasiparticle (QP) excitations are not well un- 
derstood |l|. In conventional superconductors, the study 
of non-equilibrium QP relaxation was used successfully 
to extract information on residual particle-particle inter- 
actions, as well as to pin down QP and phonon lifetimes 
(for a review, see Ref. Q)- The typical time- resolved ex- 
periment is a measurement of the change in the system's 
dielectric constant as a function of time following a pump 
pulse which creates a non-equilibrium QP distribution. 
The excited QPs decay to equilibrium over a series of 
timescales involving several steps, including at least: (i) 
a cascade of pair production until a quasi-equilibrium is 
reached between "hot" QPs of roughly the gap energy 
Ao and phonons of energy 2Ao, and (ii) slow recombina- 
tion of QPs into Cooper pairs. The timescales involved 
in step (i) are O(ps), but can be much longer in step (ii) 
[0(ns-/is)] since energy is continually exchanged between 
the electron and phonon systems until heat is removed 
at the sample surfaces; this long decay is sometimes re- 
ferred to as the "phonon bottleneck" . These two rates 
were clearly identified in, for example, recent pulsed syn- 
chrotron measurements on Pb 

A priori, several differences are to be expected in the 
cuprate superconductors. The very strong interactions in 
the normal state and the larger gap scale Ao suggest that 
electron-electron rather than electron-phonon scattering 
is the dominant relaxation mechanism. Furthermore, the 
gap anisotropy and the equilibration with thermal QPs 
at the nodes leads to a new type of slow decay, which 
we call the "antinodal bottleneck" . This is the inability 
of antinodal QPs to decay rapidly by pair-breaking while 
simultaneously satisfying energy and momentum conser- 
vation, due to the large difference in velocities at the 
nodal and antinodal points of the Fermi surface. (This 



idea has also been proposed by Orenstein 

Although several recent experiments on the cuprates 
have probed the decay of non-equilibrium QP popula- 
tions in time HEQHUEaElEl, thcir metl> 

ods and results differ considerably and it is not clear 
yet whether a consistent interpretation can be extracted. 
A short relaxation process 0(1-10 ps) has indeed been 
identified 0, 0i IS @ m experiments with pump pulse 
energy much greater than the gap scale. In addition, a 
much longer component O(10 fis) component has been 
observed [5j and identified with a recombination bottle- 
neck. Recently, however, Segre et al. 0| pointed out 
that the early relaxation component in thcir measure- 
ments on ortho-II YBCO6.52 single crystals was strongly 
T and intensity dependent, and in the limit of weak in- 
tensity they fitted it to a T 3 law. At low T, this "slow" 
component was of order hundreds of ps, but still many 
orders of magnitude shorter than that seen in Ref. [f| . 

Here we present a theory of the single aspect of the 
complicated non-equilibrium physics of pump-probe ex- 
periments most peculiar to the d-wave superconductor, 
namely the mechanism intermediate between steps (i) 
and (ii) whereby hot QPs scatter through the antinodal 
bottleneck before recombination (processes on longer 
timescales have been discussed in Ref. ^H). Within a 
model where QPs are scattered by a simple local interac- 
tion, we show that at high T there is a fast relaxation due 
to Umklapp scattering, but below some crossover temper- 
ature relaxation is dominated by diffusion in momentum 
space along the Fermi surface from the antinode to the 
nodes. This idea is a momentum-space analog of real- 
space "QP traps" in conventional superconductors 
in the sense that there is an intermediate stage of relax- 
ation in which QPs diffuse to a region of lower gap. 

Antinodal bottleneck. — The (extrapolated) relaxation 
rate of QPs of energy ~An vanishes as T — ► in the 
experiment of Segre et al. [lOj, implying both that the 
Ortho II sample is extremely pure, so that impurities 
may be neglected, and that pair-breaking processes are 
forbidden. We now show that the latter effect is a simple 
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consequence of the requirements of energy and momen- 
tum conservation. We assume a <i-wave order parameter 
Afc = A cos 29 over a Fermi surface parameterized by 9. 
In the coordinate systems of Fig. da) the BCS dispersion 
relation can then be expanded for momenta k near the 
antinode as Ek — + 2k 2 Ep/Ak and for momenta p 

near the node as E p — 2^jE F p 2 L + ^oP 2 (with k F =l). 

For simplicity, consider an antinodal QP exactly on the 
Fermi surface, which for small transfered momentum q 
must scatter parallel to the Fermi surface, 4> ~ § — ^ [ see 
Fig. nja)]; the energy lost is then E\ ost « 2gA sm26>. 
The minimum energy required to break a Cooper pair, 
i.e., to create two nodal QPs with momenta p and —p—q, 
is 2g£ F sin|f - 9\, which for |f - 9\ > A /E F is larger 
than E\ ost . More generally, pair-breaking cannot occur 
for k < 1/v / 2(A /£'f) 2 cos 2 29, where k is the distance 
from the Fermi surface. The origin of this bottleneck is 
that the velocity of the antinodal QP is much less than 
nodal velocities. While the above analysis neglects Umk- 
lapp scattering, including it only modifies the picture 
quantitatively: the bottleneck for pair-breaking remains 
but its size now depends on details of the Fermi surface. 

How, then, can a hot QP relax, trapped as it is in the 
bottleneck? It could be scattered by phonons, impurities, 
various electronic collective excitations, thermal (nodal) 
QPs or other hot (antinodal) QPs. For the creation of 
phonons there is an analogous bottleneck, although scat- 
tering from thermally excited phonons is possible. We 
now investigate in detail the consequences of scattering 
from thermal QPs; contributions from other hot QPs can 
be separated by their dependence on the laser intensity 
(see below). While our calculation includes only QP- 
QP scattering, we expect qualitatively similar results for 
scattering from (non-critical) spin fluctuations, as has 
been shown to be the case for the optical conductivity 
up to a phenomenological prefactor |l5l Il6| . 

In the spirit of a Fermi liquid approach we assume 
the existence of well-defined QPs, d\,^ — Ukc\,^ — vuc-ki, 
where Uk and Vk are the usual Bogoliubov coefficients. At 
low T two-particle processes dominate, parameterized by 
the following Hamiltonian: 

H = E [^ P - P 3^2^344t4|4|4t C 1 ) 

+ V^r^r^dl^dl^ffdio- 

where i = ki,rij — UiVj + ViUj,?^ — UiUj — ViVj and 
the sum conserves crystal momentum. We assume that 
the effective interactions V are only weakly dependent 
on momentum (this is not justified if long-ranged spin- 
fluctuations are important). Hot QPs predominantly re- 
lax by scattering from thermal QPs, as expressed by the 
last two terms; energy conservation forbids leading con- 




FIG. 1: (a) Coordinate systems for the QP momenta (heavy 
line: Fermi surface, dashed line: superconducting gap). Mo- 
menta k near the antinodes are given in polar coordinates 
(kF + k, 6) and momenta p near the nodes in Cartesian co- 
ordinates (px,P||) relative to the node. The antinodal QP 
scatters from k to k + q, where q = (q, <j>) is in polar coordi- 
nates, (b) A typical Umklapp scattering process, (c) Typical 
magnitudes of the momentum transfer for scattering perpen- 
dicular and parallel to the nodal Fermi surface, (d) Schematic 
picture of diffusion of a hot antinodal QP towards the node. 

tributions from the first term, while the bottleneck dis- 
cussed above blocks pair-breaking processes from the sec- 
ond term. The BCS coherence factors ry , fy are of order 
unity for the relevant processes and so will be neglected. 

There are two types of scattering of hot QPs from 
thermal QPs, distinguished by their typical momentum 
transfer. We first discuss Umklapp scattering, which in- 
volves large transfered momenta, q ~ fcp, as illustrated 
in Fig. IHb). A hot QP undergoing Umklapp and sub- 
sequent pair-breaking scattering can reach the node in 
a small number of scattering events. However, due to 
kinematic constraints such scatterings are only possible 
if the thermal QP is some minimum angle away from the 
node and therefore has a minimum energy A\j. Accord- 
ingly, the corresponding scattering rate is exponentially 
suppressed. Using Fermi's golden rule we obtain 

J_~ M faAvT)*/* — An/T M 2 <X 
ru qq A Q /E F ' qq " h{2i:hv F y ' 

(2) 

where a is the lattice constant, v F is the Fermi velocity 
and the size of Ay < A depends on the details of the 
Fermi surface shape. An analogous Umklapp gap deter- 
mines the momentum relaxation of nodal QPs |16| and 
hence the optical conductivity at low frequency. 

At low T the second type of process becomes impor- 
tant: normal scattering, which we now show involves 
small momenta, q <C kp. For a given q the slow antinodal 
QP loses energy much less than qAo, so that the energy of 
the fast nodal QP changes little, E p w E p ^ q . Therefore 
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the typical q depends on direction [see Fig. dc)] , and is 
q w 2p|| ~ T I Ao for scattering parallel to the Fermi sur- 
face at the nodes (<f> ss j and symmetry-related angles), 
and q ~ 2p^ ~ T/E-p in perpendicular directions. 

Since q <C fcp, the hot QP experiences many colli- 
sions with the nodal QPs before it reaches the nodes, 
as sketched in Fig. 0d). This has two important con- 
sequences: (i) the trajectory of the hot QP is diffusive 
in momentum space, and (ii) the hot QP thermalizes in 
the radial direction with the nodal QPs (see below). The 
typical momentum of a hot QP is thus determined by 
E k - A cos 29 ^ T and so k ~ k th = y/TA^/E F . 

Boltzmann equation. — To calculate the time tq for a 
dilute distribution gk(t) of antinodal QPs to reach the 
nodes we first consider the situation where the radial 
distribution of hot QPs lies entirely within the bottle- 
neck, fc th < (A /E F ) 2 , which requires T < A^/Ep. We 
assume the nodal QPs are always in local thermal equi- 
librium, as Tg 1 will turn out to be much slower than the 
nodal relaxation rate ~M qq T 3 E F / Aq. Linearizing in gk 
the Boltzmann equation reads 

dgk 



dt 



dq 



gk-qGqq(q, £k-q,q) ~ SfcGq q (<7, 



G qq {q,e) = M qq J dp/ p (l - f p - q )5(e p - q + e) 



(3) 
(4) 



Here f p is a Fermi function, Ek.q — Ek — Ek+ q is the 
energy transfer, and G qq is the integral over momenta 
near all four nodes of QP-QP scattering processes. 

Since we expect the hot QPs to thcrmalize in the radial 
direction we make the transformation 
C k (t) 



9k{t) 



-fiE h 



(l + fc)dfce 



-PE h 



(5) 



where Ckit) is a smooth function of k to be determined. 
[The factor of (1 + k) comes from the use of polar coor- 
dinates.] Using detailed balance and (JSJ), the Boltzmann 
equation now takes the form 



1 dC k 



he^ dt 



J dqG(q,e k , q ) 



C, 



k+q 



Cfe 

he,. 



(6) 



If the hot QPs are close to radial equilibrium we ex- 
pect Cfc to depend only weakly on k and therefore ex- 
pand it as C k (t) = C$ k (t) + {k~ (k)) C\ k (t) H , where 

(k) = VfcJX 1 + k)&kke~P Ek such that / dk g k (t) = 
J d9Cg(t). The equation for Cg can then be projected 
out from (0) by integrating over radial directions. Since 
qde(Cg /he) <C C'e/he we can perform a gradient expan- 
sion, leading to (for 9 <C \ ) 

d 
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\c e(ty 




dt 


he 
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(7) 



1 



-0E k 



F=- I d 2 qdk'— f — qsm(9 + q>) G qq (q, e k , q ) , (8) 



D = / d 2 qdk 



e -0E k g 2 sm 2^ + ( 
~hg 2(1 + fc) 



G m (q,s k<q ). (9) 



The main contribution is from <fi « j leading to "zig-zag" 
diffusion in momentum space [see Fig.^d)], and hence 

7^ M qq T 4 E F n=p T 
80 



F 



A 3 



4A 



(10) 



We have dropped terms in Cg (i > 1) on the right-hand 
side of since they decay faster than Cg (see below). 
As F/D = 4A /T we can rewrite for 9 < f as 



d t C° = -F d e (9C°g) + D d 2 C° e 



(11) 



Now D can be identified as the diffusion constant in mo- 
mentum space, while F parameterizes the force driving 
hot QPs towards the nodes. We emphasize that the forms 
of and (fTTJ) are completely determined by detailed 
balance, local conservation of the number of hot QPs and 
symmetry under 6 — > —9. The Boltzmann distribution 
Cg oc hg obviously solves (0), while the conservation of 
the number of hot QPs, Ah ot = J d9Cg(t), is manifest 
in fTT|) : the equivalence of and ifTTJl fixes F/D. 

For an initial Gaussian distribution of width a the 
solution to (|llfl is a Gaussian of width cr(t), where 
a 2 {t) = (2D/F + a 2 )e 2Ft - 2D/F. For short times, 
t <C 1/F, the motion is diffusive, a = 2^/Dt, while for 
t ^> 1/F we obtain a oc e Ft . The antinodal QPs reach 
the node after a time t such that a(re) ~ 7r/4, and hence 

I/to ~ F/lnimml^/F/D^- 1 }) ~ M qq T A E 2 /A 3 . (12) 

Details of the initial distribution only give log corrections. 

To check the assumption of radial equilibrium we cal- 
culate perturbatively the evolution of the corrections 
Cg(t),Cg(t) . . . using ©. Deviations from radial equi- 
librium decay rapidly, dCg/dt « ~G\/ri where r, < 
TeAg/Ep <C Te for i > 1, and can therefore be neglected. 

Higher temperatures. — We now consider higher tem- 
peratures, A^/E F <C T <C Aq/^f, where the upper limit 
comes from the neglect of curvature in the nodal disper- 
sion E p . In this regime the typical antinodal velocity 
4k t hE F /Aq = AEf ^/T/A is now large compared to the 
nodal velocity parallel to the Fermi surface, but still small 
compared to the perpendicular nodal velocity. This has 
two consequences, (i) Pair-breaking (and recombination) 
is now allowed, as k t h ^> Aq/E f , i.e. the QP is outside 
the bottleneck; in such a process one QP turns into three, 
of which two must reside at the nodes (due to energy 
and momentum conservation), (ii) Typical transfered 
momenta still turn to be small, for both pair-breaking 
and QP-QP processes, although now q ~ y/TAo / E-p for 
scattering parallel to the nodal Fermi surface (instead of 
q ~ T/A when T < Ag/S|). 

All physical requirements which led to the form of 0, 
(|llfl and fixed the ratio F/D are therefore fulfilled: con- 
servation of antinodal QPs, small momentum transfer 
and detailed balance (since both pair-breaking and re- 
combination are possible) . Generalizing (JSJ , we find 



1/T fl w (0.92A/ qq + 0.67M pb ) JT^A 3 /E F , (13) 
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with M p b defined analogously to M qq , Eq. J2J. 

Discussion. — We have shown that for vanishing laser 
intensity there are two relaxation mechanisms for hot 
QPs: Umklapp scattering from thermal QPs with an 
exponential relaxation rate ©, and normal scattering 
with a power-law rate (|12II13|I . While early experiments 
were fitted to a T 3 law [lfj, later results 0, [TlJ sug- 
gest activated behavior for the lowest measured tem- 



perature (down to 10 K), t 



je E ^ t / T , with 7 



8x 10 s , -E act ~ 100 K. This is consistent with the pic- 
ture given above, where from (|2J) and (|13|) the crossover 
to T 5 / 2 behavior is expected to occur when q- e ^/ t ^ 
TAo/£|, or T ss 8K, assuming A w 200 K and 
E-p w 8A The surprisingly low crossover scale 

^0.1i? ac t arises since many small-momentum scattering 
events are required for the hot QP to reach the nodes 
and therefore (| 1 3D is small. Note that the crossover from 
T 5 / 2 (JT3J to T 4 (H2) is expected around T ~ A 3 ,/^ 2 , w 
IK. 

As a consistency check of our approach we now con- 
sider the absolute magnitude of the relaxation rate, for 
which we require M qq in (J2J. Assuming that a similar 
matrix element appears in the optical conductivity, we 
can use the fit of Walker and Smith to YBCO6.99 
to find M qq w 2 x 10 17 eV _3 s -1 . This gives t v w 1 ns 
at T = 20 K, to be compared to the experimental value 
of 0.03 ns. Given the possible momentum-dependence of 
Mqq, Fermi surface shape effects, the difference in doping 
in the two experiments and especially the uncertainty in 
Eact, this is not unreasonable agreement. 

Within our model we can also estimate the inten- 
sity dependence of the scattering rate, which experi- 
mentally 0, is given by r" 1 w cp almost inde- 
pendent of T, where p is the energy density deposited 
by the laser and c w 1.2 x 10 6 s _1 J _1 m 3 . Assum- 
ing a uniform distribution of hot QPs on the Fermi 
surface, we use Fermi's golden rule to find l/r n f=a 
87rM qq n h ot(^F) 4 ln(£ ; F/Ao)Ao/S|, where n hot is the 
(planar) density of hot QPs. With the above numbers, 
a comparison to the experiment would suggest that only 
~0.1% of the deposited energy is used to create nhot hot 
QPs each of energy ~A ; equivalently, assuming 100% 
energy conversion, our predicted rates are a factor 10 3 
too large. While a sizable fraction of the laser energy 
probably goes into the phonon system, we think this is 
unlikely to explain the discrepancy completely. 

Conclusions. — Using BCS theory for d-wave supercon- 
ductors we have calculated the electronic contribution 
to the relaxation of hot QPs. An analytical solution to 
the QP Boltzmann equation shows that small momentum 
scattering leads to a diffusive dynamics of antinodal QPs 
and associated power-law relaxation rates, unfortunately 
for T < 10 K which has not yet been investigated exper- 



imentally. At higher T the relaxation is dominated by 
Umklapp scattering with large momentum transfer lead- 
ing to exponential T dependence, consistent with exper- 
iments on ultra-clean samples. However, it is currently 
unclear whether our estimates for magnitudes of relax- 
ation times are correct, and improvement will require a 
more complete understanding of the relation between the 
incident photon intensity and the number of hot QPs. 

There are still many details and apparent discrepan- 
cies between the various non-equilibrium measurements 
on the cuprates which we cannot understand at present, 
and it is clear that more realistic physics must be incor- 
porated to describe the evolution of the QP distribution 
over all timescales of interest, at different dopings, and 
in dirtier samples, as well as the dynamics of the order 
parameter. We believe nevertheless that the theoretical 
analysis of these experiments is worth pursuing since they 
potentially provide a window on the frequency and mo- 
mentum dependence of the interactions responsible for 
superconductivity in these systems. 
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